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We study both theoretically and experimentally the nonlinear interaction be- 
tween an intense surface acoustic wave and a two-dimensional electron plasma in 
semiconductor-piezocrystal hybrid structures. The experiments on hybrid systems 
exhibit strongly nonlinear acousto-electric effects. The plasma turns into moving 
electron stripes, the acousto-electric current reaches its maximum, and the sound 
absorption strongly decreases. To describe the nonlinear phenomena, we develop a 
coupled-amplitude method for a two-dimensional system in the strongly nonlinear 
regime of interaction. At low electron densities the absorption coefficient decreases 
with increasing sound intensity, whereas at high electron density the absorption 
coefficient is not a monotonous function of the sound intensity. High-harmonic gen- 
eration coefficients as a function of the sound intensity have a nontrivial behavior. 
Theory and experiment are found to be in a good agreement. 
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Introduction 



The interaction between surface acoustic waves (SAW's) and mobile carriers in 
quantum wells is an important method to study the dynamic properties of two- 
dimensional (2D) systems. The SAW can trap carriers and induce acoustic charge 
transport as has been investigated in a number of systems in view of possible de- 
vice applications JL]. Also, the SAW-method was applied to study the quantum 
Hall effects || [|, |4j], electron transport through a quantum-point contact lateral 
nanostructures pfl, and commensurability effects in a 2D system J7|. However, all 
those experiments have been done in the regime of small signals and linear interac- 
tion. A recent Letter by Rotter et al. || reports strongly nonlinear acousto-electric 
effects in a 2D electron gas (2DEG), which become possible in hybrid structures 
based on semiconductors and LiNbO^ [9-11]. In these experiments an in- 

tense SAW breaks a 2DEG into moving electron stripes and all characteristics of 
the acousto-electric interaction are strongly modified as compared to the linear case. 



In modern hybrid structures [|IT] the SAW-induced potential amplitude can become 
comparable with the band-gap of a semiconductor. The previous Letter || on 
nonlinear effects in the hybrid structures with a 2DEG includes a brief qualitative 
analysis. Here we present a detailed theoretical study of nonlinear acousto-electric 
effects in a 2D electron system and develop a coupled-amplitude method for intense 
SAWs interacting with a 2DEG. Using our theoretical results we can explain main 
experimental observations. Including the effect of electron diffusion we find a good 
quantitative agreement between theory and experiment for the case of the SAW 
absorption coefficient in the nonlinear regime. 

Nonlinear acoustic waves in bulk piezocrystals with free carriers were discussed in 
a number of theoretical papers [12-16]. The application of a dc voltage to the crystal 
can result in the current amplification of sound |12], |lj]] and in the formation of 
stationary nonlinear waves [|12|, 14, 16|. Analytic results in the theory of nonlinear 
acoustic waves in bulk piezocrystals with free carriers were mostly obtained in the 
limit of small amplitudes or for the case of very intense acoustic waves [113, [TBI ITB 
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For nonlinear SAW's in crystals with a 3D electron gas, a theory was developed 
in the limit of very high amplitudes, when the SAW bunches electrons near the 
crystal surface [I8|. Another theoretical aspect related to the generation of the 



second harmonic of a SAW was studied in Ref. [19] by using the coupled-amplitude 
method and perturbation theory A theory of acousto-electric interactions in a 2D 



electron system was developed mostly for the linear regime of interaction [EG, 21 



Here, we study both theoretically and experimentally the transition from the lin- 
ear regime of the acousto-electric interaction to the limit of strongly nonlinear effects 
in a 2DEG. Our theoretical results are applied for a description of experimental data 



on hybrid structures || [TTJ] . We pay attention to density-dependences of the absorp- 
tion coefficient and the SAW-velocity shift because the electron density is a tunable 
parameter in experiments on 2D systems. Such dependences were not discussed in 
detail in the context of 3D systems [12-16]. It turns out that for low densities, the 
absorption coefficient is a decreasing function of the sound intensity caused by the 
trapping of electrons in the SAW piezoelectric potential. At sufficiently high elec- 
tron density, however, the absorption coefficient is a non-monotonous function of the 
sound intensity. This behavior can be understood in terms of a dynamical screen- 
ing effect. Also, our quantitative analysis shows that the absorption coefficient at 
room temperature is strongly reduced due to electron diffusion. A non-monotonous 
behavior with increasing sound intensity was also found for the intensities of higher 
harmonics in a short device. 

The coupled-amplitude method with the introduction of fast and slow variables 



was used before for a description of bulk systems |T^, However, the formulas 
from the bulk theory can not be directly applied to surface waves because of the com- 
plex character of the lattice vibrations in a SAW. To develop the coupled-amplitude 
method in a 2D system, we introduce a local velocity shift and a local absorption 
coefficient in an integral form through an electron current, a SAW potential, and an 
electro-mechanical coupling coefficient of a microstructure. The resulting formulas 
can then be applied to any type of SAW's interacting with a 2DEG. Moreover, using 
our approach, we can find solutions for higher harmonics for arbitrary SAW inten- 
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sity. This is in contrast with perturbation methods developed earlier for 3D systems 
and valid at small SAW intensities JT^]. At very high intensities, we find analytical 
asymptotic dependences for the high-harmonic intensities in a 2D system. 

As has been shown before, the linear approximation holds when 5n <C N s , where 
N s is the equilibrium 2DEG density and 5n is a density perturbation due to a SAW. 
In the linear theory, the absorption coefficient T and the SAW-velocity shift due to 
the 2DEG 8v® are given by the well-known relation [pO , 



Sv° a J = K 2 eff (q) ia Q /g m 

v° s 2q 2 l + i(<T /a m + D e q/v°y [ ' 

where q is the SAW wave vector, is the sound velocity in the absence of a 2DEG, <xo 
and D e are the 2DEG conductivity and the diffusion coefficient, respectively. K^ff is 
the effective electro-mechanical coupling coefficient. cr m (g) = t>°e e //(g)/(27r), where 
e e ff(q) is the effective dielectric constant in a 2D system. In most piezoelectric 
crystals, the coupling -C 1. In GaAs, K^j- = 0.00064, whereas in the hybrid 
structures as studied here, it is two orders of magnitude larger, in the range of 



0.01 — 0.05 [10], but still much less then unity. The goal of this paper is to describe 



the acousto-electric interaction for the case of large amplitude SAWs, when 5n ~ N s 
and the perturbation theory is not longer valid. At the same time, the coupling K^ff 
will be assumed to be much less than unity. Below, we will generalize the results 
following from Eq. [I] for the strongly nonlinear case, when 5n ~ N s . 

The paper is organized in the following way. In the first section, we will give the 
general equations for SAWs on a piezoelectric crystal. The second section is devoted 
to a coupled-amplitude method developed for the case of intense SAWs. The third 
and fourth sections are about phenomena related to large-amplitude SAWs in a 
2DEG. Then, we will discuss experimental data on the hybrid structures and apply 
our theoretical results for the interpretation of experiments. 
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1 Model and general equations 

In usual GaAs-based microstructures it is very difficult to realize SAW's with high- 
amplitude potentials because the electro-mechanical coupling in GaAs is relatively 
weak. A strong piezoelectric interaction can be achieved in hybrid structures || [9| 
(Fig. 1). Those consist of a semiconductor layer being bonded to a piezoelectric host 
crystal, in our case, LiNbOs || |9|] . The semiconductor layer contains an InGaAs — 
AlGaAs quantum well (QW) with a high-quality 2DEG, to which Ohmic contacts 
are formed. The distance between the QW and the piezocrystal is only 32 nm, 
whereas the distance between the QW and the top transport gate is d = 450 nm. In 
our model x and y are the in-plane coordinates and z is the normal one. The QW 
plane corresponds to z — d, and the SAW travels in the x-direction (see Fig. 1). By 
changing the transport-gate voltage Vt one can tune the electron density N s (Vt) in 
a QW. In this structure, traveling SAWs can induce very strong piezoelectric fields 
in the semiconductor layer due to the strong piezoelectricity of the host LiNb0 3 
crystal. A SAW is induced and detected by the metallic interdigital transducers, 
I DTI and IDT2, respectively, at room temperature ||. The acousto-electric current 
is measured between two Ohmic contacts labeled 1 and 2 in Fig. 1. 

SAWs in a piezoelectric crystal with a 2D plasma are described by the system 
of equations |22| 



put = Cikimd m d k ui + pukdid k (j) (2) 
e(x 3 )didi(p - A-npikididiUk = -4nen5(x 3 - d), (3) 

where e = — |e| is the electron charge, p is the mass density, and e is the dielectric 
constant. Further, Cikim is the elastic tensor, pi ik is the piezoelectric tensor, u k (x, z, t) 
is the lattice displacement, and <p(x,z,t) and n(x,t) are the electrostatic potential 
and the 2D electron density, respectively. In Eqs. |],[| we have used the notations 
dif = df /dxi and / = df /dt, and the sum convention for repeated indexes, t is the 
time, and x 3 = z, x 2 = y, and X\ = x are the coordinates. 

In our geometry, the SAW propagating in the x-direction is a purely Rayleigh 
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wave, in which only two components of the displacement, u x and u z , are nonzero. 
Hence, the electric field E is also polarized in the (x?/)-plane. This case corresponds 
to the hybrid structures studied in experiments [§, [TO] , where the 128°-rotated 
y-cut of LiNbOs is used. The surface of the thin GaAs film is (001). The SAW 
propagates in the [110]-direction of GaAs and X-direction of LiNbO^. 

The electron 2D plasma is described by the usual hydrodynamic equations 



en(x,t) + ^^ = 0, (4) 

j = _ eP^fo'J, (5) 

ox ox 

where j is the 2D electron current, a = \e\fm s (x,t) is the 2DEG conductivity, \x is 
the mobility, and D e is the electron diffusion coefficient. Eq. ||] is valid in the long 
wavelength limit, when ql e <C 1, where l e is the electron mean free path. 

The wave equations (|2|,|3|) should be solved together with the boundary conditions 
at the surface z = and at the interface z = d. At z — 0, these conditions are the 
following: (f)(z = 0) = and a zi = c zUm (d m ui + diu m )/2 + pi zk d t (f) = 0. Here a zi is 
the z-component of the stress. The top metal gate is thin and does not influence 
the boundary condition for the stress tensor. At the semiconductor-piezocrystal 
interface, 0, u iy and a zi should be continuous, and D z [z — d + 5) — D z (z — d — 
5) = 4iien, where 5-^0 and Di = —e(z)di(p + 2irp ik i(dkU[ + diUk) is the electric 
displacement. For simplicity, we assumed above that the 2DEG is located directly 
at the semiconductor-piezocrystal interface. 



2 Coupled-amplitude method 



The system of nonlinear equations (@-|>D can be simplified in the limit of weak electro- 
mechanical coupling, K^ff ~ p 2 /(ce) 1. In this limit, we can introduce two 
coordinates, the "slow" variable x and the "fast" variable X\ = x — v° s t |TJ, [T5| . The 
solution is of the form (see Appendixes 1 and 2), 



6 



u(x, z,t) = u(x,xi,z) = u (x, z) + a n (x)\]°[z] q n + 5q n (x)}e ignX1 + c.c. (6) 

71=1,2,... 

4>(x,z,t) = (f>(x,x 1 ,z) = (f> (x,z)+ a n (x)<S>°[z;q n + 5q n (x)}e iqnXl + c.c, (7) 

n=l,2,.. 

This solution is written as a sum of harmonics with wave vectors q n = nq, where 
q > is the wave vector of the initially generated SAW near x = 0. The vec- 
tor A°[z;g] = (u [2:; q}; §°[z; q]j and the quantity 5q n are determined by a linear 
system of equations as given in the Appendix 2. The envelope functions a n (x) 
are slowly changing on the scale of A = 2n/q. It is assumed that the SAW in- 
tensity related to the vector A°[,2;g] is unity and thus the total SAW intensity is 
haw — J2 n =i,2,..\ a n( x )\ 2 - The functions u (x, z) and (p (x,z) describe the static 
spatial distributions, that can be induced by a SAW. 

The functions j(x,xi,t) and n(x, x±,t) can be written in the standard form 
f(x,xi,t) = fo(x) + J2n=i,2.... fn(x)e iqnXl + c.c, where f n (x) is an envelope func- 
tion. Also, the electric field E can be written in the way similar to Eq. [7[ 

The solution (||[7|) is a sum of linear-like SAWs, which slowly vary in the am- 
plitude and in the ^-profile. The ^-distribution of lattice displacement U° is a sum 
of exponential functions exp [— jj(q n + 6q n (x))z], where 7j(g) are the coefficients de- 
pending also on material constants |P2fl . On short distances (~ A) the envelope 
functions a n (x) can be regarded as constants and we can solve the equations (|||-[5D 
considering only the "fast" variable x\. n n (x) and <fi n (x) should be considered here 
as the functions of the parameters a 1? 02, Uq, and E . Then, having n n and <p n as 
functions of a\,a 2 , n , and E Q , we can find the behavior of a n (x) on long-range 
scale, x ~ l/r° > A. 

The electrostatic potential is written in a self-consistent way: 

<f>(x,x 1 ,z) = (f> ind + ( f> saw , (8) 

where <p ind and <p saw are the potentials induced by a 2DEG and by piezoelectric 
charges of a SAW. Using Eq.^| we write 
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e{z)d i d i (j) ind = -AnenS(z - d) (9) 
e{z)d l d4 saw = 4-KPikidAuk. (10) 



(ft and 4> saw can be expressed by the harmonic amplitudes and 0* au; . For exam- 
ple, <p ind = <j)™ d (x, z)+E n =i A ... <Pn d (%, z)e ignXl +c.c. In the limit r max A ~ K 2 eff < 1, 
we find from Poisson's equation (see Appendix 1) and from the conservation of charge 



C d (x, d) = 27r 7 n f } , E Z d & d) = -iq n ^ d {x, d), (11) 

e eff{Qn)<ln 

j n (x) = v° s en n (x), 

where n = 1,2,3... . Here e e //(<?) = [e P + e s coth(\q\d)]/2 is the effective dielectric 
constant including the gate electrode effect, and e p and e s are the dielectric constants 
of a host piezocrystal and a semiconductor, respectively |23| . 



The nth-harmonic of the SAW-potential cj) s ^ LW is given only by a n (x) and by 
material constants and can be easily found from the Poisson equation (|I0D. At 
z = d we have 0* au, (a;, d) = C n a n (x), where the coefficient C n depends on the 
geometry and the material constants. For example, in a crystal of the type of 
GaAs, C n = P4,g(q n ), where g(q n ) is a complicated function of q n . Below, we will 
give the necessary relations for the hybrid structures. 

To find the harmonic amplitudes n n and n , we have to solve Eqs. |3|, f|, and 
|5] regarding the "slow" variable x as the constant. The slowly- varying quantities 
n n (x) and j n (x) can be found as Fourier components of the solution n(xi,x) from 
a "fast" equation in terms of X\. In a self-consistent approach, the electron density 
n(xi,x) is determined by the SAW-induced potential at z — d, <j) saw (xi,x,d) = 
4>Q aw (x,d) + J2n=i,2,... C n a n (x)e iqnXl + c.c, through a non-linear equation ||. By 
using Poisson's equation, the results of Appendix 1, and Eqs. [| and ||, we obtain 



\e\n{x 1 ,x)fx[-J-{l dx'^xr-x'^n^x^ + E^^x)] (12) 



dx\ 



dn{x u x) 

-eD e ev s0 n{x 1) x) = b , 

dx\ 
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/+oo ^—ik(xi—x' 1 ) 
dk— -—, (13) 
\k\e ef f(\k\) 

where 60 is a constant, which occurs after one integration in the conservation-of- 
charge equation (Q). Taking into account r°A < 1, we assume that the solu- 
tion of Eq. [12] is periodic in x\, n s (xi,x) = n s (xi + X,x). Also, < n s (xi,x) >= 
Jq n s (xx, x)dxi/X = n (x), where n (x) plays the role of a "local" 2D density. 

The functions a n (x) are connected by the system of nonlinear equations (see 
Appendix 2) 



^ = iSq n [a 1 , a 2 , ...;n , E )a n . (14) 
ax 



where 



r r K lfMn) 2tt en n {x) . Kl ff {q n ) ]n {x)E s n aw * 
dq n \ai, a 2 , n , E Q \ = — ; — — = 1 J - J — — ; - - 15 

and n = 1, 2, 3... . 

The local velocity change of a SAW 5v n and the local absorption coefficient T n 
can be expressed by Sq n : 



oq n = -q n — + i—- (16) 



Using the relation 



K 2 eff (q n ) = 

(see Ref. [23]) and Eqs. [TI] we write 



2\E, 



saw 1 2 



Om(<?n) 



qnln 



(17) 



fo n (x) _ < ^ w {x 1 ,x)j{x 1 , x) > ^ 



V 



Here we use the notations < f(xi,x) >= J A /(xi, x)dxi/X, and f n (xi,x) = 
f n (x)exp(iq n x 1 ) + c.c. . 
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In this section, we have assumed that ^f// ^ 1 anc ^ neglected the terms 
d 2 a n /dx 2 oc 5q n da n /dx oc -K'e// an d d5q n /dx oc -K'e//- The static electric fields 
4>o(x, z) in Eq. [7], which can be induced by a SAW, will not play an important 
role in this paper because we will consider the case with no voltage applied to the 
Ohmic contacts (Vi = V2) and a relatively short device with L <C 1/T . Thus, 
<p (x,d) ~ const and the static electric fields E 0x (x,d) can be neglected. Eqs. [14] 
and [15] can now be applied to various types of SAW's by introducing specific electro- 
mechanical coupling coefficients. 



3 Acousto-electric transport in a two-dimensional 
plasma 

Here, we intend to consider the SAW absorption and the acoustic charge transport in 
the regime V\ — V2. The behavior of a 2DEG in an intense SAW at small distances 
(~ A) can be assumed to be periodic and is described by nonlinear equation (|12|). At 
long distances (^> A) the SAW behavior is determined by complex amplitudes a n (x), 
that can be found by Eqs. [TJ| and |T5|. First, we define the boundary conditions at 
x = 0: oi(0) = y/h{0) and a„(0) = for n = 2, 3, .... Here, I x (0) denotes the SAW 
intensity generated by IDT1. In this case, for a relatively short sample the SAW 
contains mostly the fundamental SAW harmonic n = 1 with a small admixture of 



the higher harmonics n — 2,3, .... So, it follows from Eqs. [L4 and 15 



ai ( x ) = ai(0)(l +ix5q 1 [a 1 (0),0,0, ...;no,0]). (20) 

We assume that [/i(0) — h(L)] <^ /i(0), where L is the length of a semiconduc- 
tor film and I\{L) is the intensity detected at the IDT2. The latter is valid in 
a short sample where T\L 1. The first-harmonic absorption coefficient per 
unit length, which is measured in the experiments, is [ii(0) — ii(L)]/(ii(0)L) = 
Im(5q n [ai(0), 0, 0, no, 0]) = T 1 (a 1 , 0, 0, n , 0) = r x and the velocity shift is 
5vi(ai, 0, 0, n , 0)/v° = Re(5qi[ai, 0, 0, n , 0]/q) = 5v 1 /v° s . 
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The amplitudes of higher harmonics in a short sample turn out to be given by: 



a n (x) = ix \hna n 5q n [a 1 (0),0,0,..,a n ,...;n ,0]=ix e //( g ") en ^\ (21) 

a «^° € eff{(ln) 

where n = 2,3,4,... . Here, we assume that [ii(0) — h(L)] ^> I n , which is valid 
for a short sample where TiL <C 1, and we take into account that for the strongly 
nonlinear case Ti ~ T n . The intensity of the nth-harmonic is given by I n (x) = 
\a n (x)\ 2 and 

lr |2 _ \<Pn W \ 2 Kl ff {q n ) 
\a n \ 2 2q n a m {q n ) 

(see Eq. |T7|). By using Eq. ^T] we find the intensities of high harmonics at x = L: 
I n (L) = lim an _ \aJq n [ai, 0, 0, a n , 0, 0] \ 2 L 2 = nL 2 (K 2 eff /t eff )q n v Q s e 2 \n n {$)\ 2 . 

To calculate Ti and 8v\, we numerically solve Eq. [12] for the parameters a\ ^ 
and a n = when n — 2, 3, .... Then we can find the Fourier components of n(xi, a\). 
We now calculate the electron density n(x\, a\) in the long wave-length limit qd -C 1. 
In this limit the function given by Eq. [1^ is reduced to G(x\ — x[) = (4:ired/e s )5(xi — 
x[) and Eq. [12] now writes as 

... . , A-neddnixi, aA ,_,„„,.. . dn(xi,ai) 
\e\n(x u ai)u{ y ^ 11 + F^ aw cos{qx 1 + Xi)} - eD e v > (23) 

— ev °n(xi, ax) — b = 0, 

where the SAW-induced piezoelectric field is taken in the form E saw = F[ aw cos(qxi + 
Xi), and Xi is the phase of the first SAW harmonic. For simplicity, we put in the 
following Xi = 0. The constant b in Eq. |23] is directly connected with the kinetic 
motion of a SAW and thus vanishes for i>° — ► 0. Dividing Eq. |23] by n and then 
integrating over x± we get 



Aired Ff aw 
|e|u{ n(xi) H sin (qxi)} — eDJn[n(xi)] = (24) 



r x i dx\ 



consti + ev° s x\ + bo 

Jo 



i 



n(xi 
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The left-hand side of the above equation is periodic and thus the right-hand side 
should be periodic as well. It implies that 60 J A = — ev®\. 

A numerical solution of Eq. |23| n(xi) at T = 300f^ and at various averaged 
densities N s =< n(xi) > is shown in the inset of Fig. 3. Here, we use the following 
parameters: A = 33 fim, fi = 5000 cm 2 /Vs, D = fi(KT/e), e s = 12.5, and v® = 
3.8 * 10 5 cm/s. It is seen that that with decreasing N s the formerly homogenous 
2DEG turns into moving electron stripes. 



It follows from Eqs. [18,19 



8v\ < &l aw sin(qxi)j(xi) > < F( aw cos(qxi)j(xi) > 

If = 2^ ' 1 = h ' ( } 

The calculated absorption coefficient I\ as a function of the electron density N s for 
various potential amplitudes ^ s l aw = F( aw /q is shown in Fig. 2. It was calculated 
from Eqs. ^5| and the numerical solution for n{x\). The potential amplitude $* a ™ 
can be easily connected with the SAW intensity I\ by using Eqs. [T7| and E2I with 
Keff found numerically in Ref. 0. As an example, in the hybrid structures at 
SAW frequencies / = 114 MHz and 340 MHz, K 2 eff is about 0.015 and 0.035, 
respectively. 

We see in Fig. 2 that with increasing $^ aw the absorption coefficient in general 
decreases and its maximum is shifted to the higher values of N s . This nonlinear be- 
havior can be understood qualitatively as follows. For the densities N s iV™ aa: , the 
electron plasma forms the moving charge stripes and the electron velocity j/ (eN s ) 
is very close to its maximum i>°. Here i\T™ ax denotes the density corresponding to 
the maximum of the function Ti(N s ) at fixed $5""". We will show below that in the 
case N B < N™ ax Ti ~ \e\N s (v®) 2 / (hfi) oc N s /Ii . This asymptotic behavior follows 
from the Weinreich relation |p5| . In the region N s 3> N™ ax , the electric current aE x 
decreases with increasing N s because of the screening effect. Thus, the absorption 
coefficient decreases as well. 

In Fig. 3 we plot T 1 as a function of the SAW potential &\ aw for a fixed 2DEG 
density. For electron densities less than about 10 10 cm _1 , the function Ti(§l aw ) is 
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always decreasing. At higher densities Ti(^l aw ) has a maximum. We attribute this 
behavior to the screening effect in a 2DEG. At high density and small $f a ™, the 
absorption is strongly suppressed because of screening. The intense SAW, however, 
modulates the 2DEG and, consequently, reduces screening. Thus, the absorption 
coefficient starts to increase with increasing $^ aw , when $* aui is not so large. At 
larger $>* a ' w j the plasma becomes broken into stripes and T 1 decreases with Ii as 
1/h li. 

The velocity shift 5v\ numerically calculated by Eqs. |25| is shown in Fig. 4. 
With increasing SAW intensity 5v i increases and approaches to the velocity v® that 
corresponds to the case of a totally depleted 2DEG. By analyzing Eqs. ^5] and the 
function j(x\), we can see that, in the limit I\ — > oo, 5v± oc l/$* aw oc 1/v^T, like in 
a 3D plasma |T2] ]. 

In Fig. 5 we show results of numerical calculations for the intensities of the high 
harmonics I n /I n , ma x = \n n \ 2 /N^, where 

In, max = ^ Oo) = TtL 2 ^q^N* . 

€ eff 

This formula was obtained taking into account that \n n \ — > N s when $* aw — > oo. 
The quantities I n were calculated from the Fourier components of the function n(xi). 
At small amplitudes $* a ^, I n oc (& s l aw y n , being typical for weak nonlinearity. At 
larger c^ 1 ™, the behavior of I n is quite complex and strongly differs from the weak- 
nonlinearity behavior. At very large $* aM, ; the high-harmonics intensities tends to 
saturate: I n — > I n ,max- At higher electron densities the saturation of I n occurs at 
larger SAW potentials (see Fig. 5). 

4 Analytic results for acousto-electric effects in 
a 2D electron plasma 

In this section, we will give some analytic expressions describing various acousto- 
electric effects in a system with a 2DEG. For simplicity, we will not take into account 
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diffusion, hence assuming D e = in the formula for the current fl5|). This assumption 
is well justified at large electron densities or at low temperatures. It is convenient to 
start with the Weinreich relation for a 2D system in the nonlinear regime JT2], Z5 |. 



Assuming that n{x\) and j(x\) are periodic functions we can rewrite Eq. 12 in the 
form 



j(xi) -v° s en(xi) =b . (26) 

Obviously, it is valid in the limit i^f// ^ 1- Using Eq. |26| we can write < j >= 
\e\fi < nE x >= —fi < jE x > /v° s . From this equation we now get the Weinreich 
relation: 

<JE*> = _vj 
<]> V 

where < jE x > is the dissipation in a SAW. To get Eq. ^7], we have neglected the 
averaged electric fields < E x >= assuming the case V\ = V?- 

We now consider Eq. ^3| with D e = 0. At fixed $* aw , there is a critical density, 
N^t, for the formation of stripes in a 2DEG. When N s < N cr u the plasma is split 
into electron stripes. If A^ s > N crit , the plasma is continuous but can be strongly 
modulated in space. For the case N s < N crit , the constant bo in Eq. |23| becomes zero 
and thus Eq. ^3] has formally two solutions: 

e s ^ s , aw v° e s 

fi(xx) = const 2 - , | | - sin (qx{) s x; h{xi) = 0. (28) 

47r|e|a (i 47r|e|c/ 

The solution has to be a periodic continuous combination of these two functions. 

It follows from Eq. |28| that in the limit cr ^> a m N crit ~ $f om e s /(47r|e|d), where 

(T = \e\(J,N s . 

When N s < N cr u the plasma turns into stripes, which means that electrons 
are totally trapped and the local electron velocity in a 2DEG reaches its maxi- 
mum j/(eN 8 ) = v®. From the Weinreich relation, we get r x =< jE x > jl\ = 
\e\N s (v° s ) 2 /{»h) oc l/h 1,0. 
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Eqs. ^ and |24] can be used to find an asymptotic formula for I\ and < j > 
in the large density limit when o~o ^ &m- In the limit o~o ^> <J m , n ^ N s — 
(e s ^ s 1 aw )/(An\e\d)sm(qx 1 ) (see Eq. By using the above and Eqs. and |27|, 

we have in the limit a 3> Cm and in the region $l aw < § cr u'- 

* = (29) 

where T max = and & cr it = A^" s 47r|e|<i/e s . This equation is valid when ($ C r«t — 

$f aw )/$ crit > a m /<r . Eq. |29| reproduces the numerical data for Tif^f"") in Fig. 3 
at large densities. The asymptotic formula (^) was given before in Ref. f26[ without 
noting the condition cxq ^> a m . In the linear regime of interaction I\ and 5v\ are 
given by the formula ([j]). 

In the end of this section we consider an asymptotic behavior for the high- 
harmonic intensities /„ in the limit I\ — > oo. The electron density n s {x\) at high 
Qsaw can k e wr i^ en j n a parabolic approximation. Then, by calculating the Fourier 
components n n , we find that I n ,max — I n oc ($^ aw ) -2 / 3 oc I x 1 , where n = 2, 3, .... 

5 Comparison with experimental data 

The experiments involving SAW's were performed on the hybrid semiconductor- 
LiNbOs structures fabricated by the epitaxial lift-off (ELO) technique developed 
by Yablonovich et al. |27|. The structures contain a 12 — nm-thick high-quality 
Ino.2Ga . 8 As quantum well (QW) embedded in modulation doped Al 02 Ga ,sAs bar- 
riers. In these structures, the thin semiconductor layered system including a QW 
was tightly bound to the lithium niobate host crystal by the van der Waals forces 
|8|, [|, [nj. The MBE grown quantum well structure is removed from its native GaAs 
substrate by etching an AlAs sacrificial layer below the active semiconductor sys- 
tem. The thin ELO film with a thickness of only 500nm is then transferred onto 
the host LiNbOs crystal. The parameters for this structure were described already 
in the beginning of Sec. 1. The geometry of the structure is shown in Fig. 1. For 
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further details related to the fabrication procedure of such quasi-monolithic struc- 
tures we refer the reader to Refs. [8-10]. The experiments were performed for two 
SAW frequencies / = 340 MHz and / = 114 MHz at room temperature. 

The SAW in our experiments can be strong enough to break up an initially 
homogenous 2D plasma into moving stripes. The transition to the regime of mov- 
ing electron stripes was directly observed in the experiments on acoustic charge 
transport (ACT) in samples with specially-designed injection and detection dates 
||. In these experiments the velocity of the ACT-signal first increases with the 
SAW-intensity and finally saturates at the sound velocity. The latter manifests the 
formation of stripes. Strongly nonlinear effects are also observed in the attenuation 
data. The attenuation of a SAW with / = 114 MHz for different intensities is 
plotted in the insert of Fig. 2 as a function of the transport-gate voltage, which 
determines the averaged electron density in a 2DEG. At small SAW-intensities, the 
electronic sound attenuation T as a function of the conductivity <jq is described 
by the well-known linear-theory equation ([]]) and exhibits a maximum. This lin- 
ear regime is realized in our experiments at the smallest SAW intensities of about 
—12 dBm (insert of Fig. 2). It is seen from the insert of Fig. 2 that at high SAW am- 
plitudes the attenuation is strongly suppressed and its maximum is shifted to higher 
gate bias or conductivity, respectively. The experimental data for a SAW with the 
frequency / = 340 MHz look qualitatively similar to those for / = 114 MHz and 
were given earlier in Ref. ||. 

The nonlinear regime of interaction is described by the theory given in Sees. 2, 
3, and 4. To quantitatively compare theory and experiment, we now express the 
SAW potential amplitude ^ s l aw through the input radio frequency (RF) power P. 
The SAW intensity can be written as L\ = I saw = 2(P/w)10 _/L//10 , where the width 
of the transducer w = 0.55 mm. The insertion losses {I Li) in the transducers were 



measured to be 15 dB. Then, the SAW potential can be written using Eq. 22 as 



$saw = K eff ^2Lx/{qa m ), where K 2 eff = 0.015 for / = 114MHz [Kj. In Fig. 3 we 
also show the experimentally measured absorption coefficient at the gate voltage 
V t = —7.5 V. This voltage corresponds to the maximal attenuation for the smallest 
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RF power, — 12 dBm (see insert of Fig. 2). From the linear theory we find that the 
absorption coefficient is maximal at N s = 0.9 * 10 10 1/cm 2 . One can see from Fig. 3 
that the experimentally measured function T 1 ((^ s 1 aw ) for V t = —7.5 V is in a very 
good agreement with the calculated one for N s = 0.8*10 10 cm~ 2 . Here we did not use 
any fitting parameters. This quantitative agreement becomes possible if we account 
for the diffusion coefficient. The maximal absorption coefficient as calculated from 
the linear theory is T^/J x = 7.6 cm -1 . Without diffusion this value is about T max = 
14.3 cm' 1 . Thus, the diffusion strongly suppresses the SAW absorption. 

At fixed SAW power and a sufficiently small density N s , the 2DEG is divided 
into stripes and Ti increases with increasing the gate voltage. In our theory, Ti oc N s 
in the regime of stripes, which explains the increase of the attenuation at small gate 
voltages in the insert of Fig. 2. At a sufficiently large gate voltage, the absorption 
coefficient as a function of V t starts to decrease because of the screening effect in a 
high-density 2DEG modulated by a SAW. The interplay of these two effects leads 
to the shift of the maximum of the function Ti(Vt) shown in the insert of Fig. 2. 

The change of the SAW- velocity due to the electron plasma is shown in the insert 
of Fig. 4. With increasing RF power the curves in Fig. 4 are again shifted towards 
larger electron conductivity which can be understood in terms of screening. With 
increasing SAW intensity the electron plasma is strongly modulated or even split into 
stripes and the screening of piezoelectric fields by electrons becomes not so effective. 
Thus, the shift of the SAW-velocity due to electrons decreases with increasing the 
SAW intensity. The experimentally observed shift is in qualitative agreement with 
our modeling shown in Fig. 4. However, our theory does not reproduce the character 
of Svi(Vt) in the region of small electron densities. Likely, the measurement of the 
SAW velocity is not so sensitive to a low-density electron system in comparison with 
the attenuation method. 

In Fig. 6, we show the quantity YiP/I ae as a function of the RF power to verify 
the Weinreich relation. The acousto-electric current < j >= I ae (P) was measured in 
a "short circuit geometry", where the Ohmic contacts are directly (without resistor) 
connected to the current measurement instrument. For more details on acousto- 
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electric current measurements we refer to Ref. 11]. We see from Fig. 6 that the 
ratio TiP/I ae has a weak power-dependence. Thus, our experimental data are well 
described by the Weinreich relation. Slight deviations from the Weinreich relation 
seen in Fig. 6 can come from the density-dependence of the mobility /j,(N s ), that is 
expected to be relatively weak at room temperature. The reason is that the main 
electronic scattering mechanism at high temperatures is due to acoustic phonons 
and is relatively insensitive to the 2D density. 



6 Conclusions 

The theoretical results obtained in Sec. 2, can also be applied to study dynamics of 
SAW's at large distances in a long sample, where the contribution of high harmonics 
can be very important |28|. [Z9|| . This long-distance transformation to high harmonics 
was studied experimentally for SAW's interacting with a 3D electron gas of a semi- 
conductor on a piezocrystal [p8 , 29[ . In the presence of a dc voltage applied to the 



crystal, it is possibly to expect the appearance of nonlinear waves with a stationary 
profile or with a stationary energy flow [IT3, IT3, E58|. In a wave with a stationary 



energy flow the wave shape is periodically changed in space p8| . Another scenario 
can relate to chaotic dynamics in an acousto-electric system []30| . Eq. |14] can be used 



to numerically model these phenomena in 2D electron systems at long distances. 

To conclude, we have studied strongly nonlinear acousto-electric phenomena 
caused by the interaction between a SAW and a two-dimensional electron system. 
In the experimental measurements performed on hybrid semiconductor-piezocrystal 
structures the SAW attenuation, the SAW velocity change, and the acousto-electric 
current are strongly modified in the nonlinear regime due to the formation of mov- 
ing electron stripes. By using a coupled-amplitude method we have modeled the 
decay and transformation of SAW's in the nonlinear regime. Using our theoretical 
results, we could explain our experimental findings and distinguished between dif- 
ferent regimes of the nonlinear acousto-electric interaction at large SAW intensities. 
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Appendix 1 

Here we intend to briefly discuss the electrostatics of the hybrid structure. The 
spacing between the 2DEG and the top metal gate in the fabricated structures is 
much larger than the distance from the 2DEG to the AlGaAs — LiNbO^ interface. 
Thus, to model the screening effects, we will assume that the 2DEG is located 
right on the AlGaAs — LiNbO% interface. It is convenient to solve this problem 
by a Fourier transform in terms of x and by remaining the vertical coordinate z. 
The relation between the Fourier components of the electrostatic potential W md [z; k] 
induced by 2D electrons and the 2DEG density is found from the Poisson equation 
and from the corresponding boundary conditions, 



w ^ k ^w^wf^ (30) 

where e e ff(k) = (e p + e s coth(\k\d))/2. The function G{z) = f^jS, when < z < d, 
and G{z) = e~\ k ^ z - d ^ for d < z 

Taking the electron density in the form of n s (x, xi,t) = no(x)+J2 n =i,2,... n n (x)e iqnXl - 
c.c, we can write for the induced electrostatic potential: <p md (x, xi, z, t) = 0Q nd (x, z) + 
Hn=i,2,... 0n nd ( x 5 z)e iqnXl + c.c, where q n = nq. Using Eq. we find 4>™ d (x,z) = 
2iren n (x)G(z) / ^q n e e ff(q n )^j+5(f) n (x, z), where n = 1,2,.... The correction S(f) n (x, z) ~ 
da n (x)/dx ~ 5q n ~ K^j (see Eq. [15] for Sq n ) and is small compared to W ind [z; q n ]. 
Thus, regarding n n (x) as the constants, we obtain at z = d: (j) l ™ d (x, d) ~ 2iien n (x)/[q n e, 
and E nx (x,d) ~ —iq n (j)^ d (x,d). In other words, we consider the envelope functions 
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n n (x) as constants and solve Poisson's equation in terms of the "fast" variable X\. 
Again, it is valid in the limit T°/q ~ Keff ^ !• 

As an example, we now calculate the correction 5(j) n in the linear regime of inter- 
action, when n(x, xi, t) = N s +ni(x)e iq3>1 +nl(x)e~ iqxi with 711(2;) = hie~ r ° x / 2 ~ tgSv °/ v ° : 
From the Poisson equation we find S(pi(x,z) ~ (r°/2 + iqSv s /Vg)ni(x)F(z), where 
F(z) ~ 1 for z ~ 1/q. 



Appendix 2 

In order to solve a system of nonlinear equations we will use some of results 
from a linear-response theory 00, ETJ. In a linear theory the total electrostatic 
potential and the 2D density can be written as W[x, z,t] = W[z, k]e ~ iU>t and 
n(x,t) = n k e lkx ~ luJt , respectively. Here u is the SAW frequency. It is convenient to 
introduce a quantity Hfc by means of the relation erik = — HkW[d; k), where W[d, k] 
is the Fourier component of the electrostatic potential at z — d. Eqs. ||]|| are now 
written as 



u 2 pui + c mm d m d k ui + p Kk did k W[x, z, t] = 0, 
Airp ik ididiu k + (edidi - 4irll k 5(x3 - df)W[x, z, t] = 0. 



(31) 
(32) 



Above equations should be solved together with the necessary boundary conditions 



considered in Sec. 1. Then, we rewrite Eqs. pl,32 in the form 



LhnAun — 0, (33) 

Here Lu n is a linear operator and A^ n = (u[x, z, t], W[x, z, t]j = (\J°[z; k], W°[z; kfje 1 
It follows from the boundary conditions that k = q + 5qu n pi| , where q = oj/v° s and 



qK 2 eff U q /Il q o qeeff(q) , . 

5qiin - 2 i + iynS' q ~^^' { ) 



20 



Now a solution of Eq. ^3| can be written as Au n = A°[z; q + Squ n }fo(x)e iqxi , where 
fo(x) = e lSqiinX and x\ = x — v®t. A°[z; q + 6qu n ] is a vector, that can be found from 
the matrix given by the boundary conditions |2T . 

Eq. ^ contains first and second spatial derivatives and can be written as Lu n An n = 
e iqxl [f (x)L (q)A° + fl ) (x)L l (q)A + f[;(x)L 2 (q)A } = 0, where /' = df/dx. Neglect- 
ing the second derivative fd(x), that is ~ 5qf in ~ K^ ff , we have 



Lo(g)A + i«Jg Kn Li(g)A = 0. (35) 

Now we turn to the nonlinear theory, where the equation LA = is also 
valid. The operator L is determined by the equations similar to Eqs. |3ll - |3~3l but 
with nonlinear quantity H n (x), that is defined by H n (x) = —en n (x)/(f) n (x,d). In 
the nonlinear case the vector A = Aq(x,z) + J2 n =i,2 ... a n(x)A n (x, z)e iqnXl + c.c. 
Each term in the equation LA = should be zero, and so La n (x)A n (x, z)e iqnXl ~ 
e iqnXl [a n (x)L (q n )A n + a' n (x)Li(q n )A n ] = 0. To get the latter equation, we have ne- 
glected d5q n /dx ~ Ktff an d d 2 a n /dx 2 ~ K^ff- We can solve the equation LA = 
if we choose A n (x, z) = A°[z, q n + Sq n (x)], where the vector A°[z; q] is defined above 
in the linear theory. Using Eq. B3] we get 



= iSq n (x)a n (x), (36) 

dx 

where 5q n (x) is given by the equation for Sqn n (Eq. |33|) with corrections Il g — ■> H n (x) 
and 11° — > n° n = n°. Eq. ^ is used in Sees. 1,2 to describe the acousto-electric phe- 
nomena in a 2DEG. Using Eqs. ||| and 11 and the results of Appendix 1, the denomi- 
nator in Sq n (x) (see Eq. |33|) is rewritten as 1 + H n (x)/U^ = — H n (x)4>™ w (x)/en n (x). 
Then, by using Eq. |34| and the conservation-of-charge equation j n (x) = v®en n (x), 
we obtain 



X ( )_ K hf(Q") ^ en n(%) = K 2 eff (q n ) ] n (x)E S n aW * ^ 

qAX) 2 e eff (q n )^(x) 1 2Wy s (x)\* ' 1 } 
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We use this equation in Sec. 2 (see Eq. |l^). To obtain above results, we have 
neglected the terms like a"(x) and Sq' n (x) assuming that Khf is a small parameter. 
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Figure captions 

Fig. 1. The cross section of a hybrid semiconductor-piezocrystal structure. An 
epitaxial lift-off film has a thickness 0.5 /an. The Ohmic contacts are formed to a 
2D electron gas. The transport gate with applied voltage V t is used to control the 
conductivity of the electron plasma. A high-frequency (RF) signal is applied to the 
metal interdigital transducer IDT1 in order to generate surface acoustic waves. A 
surface acoustic wave propagates though a sample and is detected by the transducer 
IDT2. 

Fig. 2. The calculated absorption coefficient of the first harmonic as a func- 
tion of the carrier density N s for various potential amplitudes $l aw . A = 33 fim, 
pi = 5000 cm 2 V/s, and T = 300 K. Insert: The measured attenuation of a SAW as 
a function of the gate voltage for different high-frequency (RF) powers applied to 
the IDT1; / = 114 MHz. 

Fig. 3. The calculated absorption coefficient of the first harmonic Ti as a function 
of the potential amplitude §l aw induced by a SAW for various fixed densities N s . 
The parameters are similar to those in Fig. 2. The dots show the experimentally 
measured absorption coefficient r\ at the gate voltage —7.5 V. This voltage corre- 
sponds to the maximal attenuation at the smallest RF power. In the inset we plot 
the calculated local carrier concentration n as a function of the in-plane coordinate 
X\ for different total carrier concentration N s . The numbers attached to the plots 
correspond to N s in units of 10 10 cmT 2 . 

Fig. 4. The calculated SAW-velocity change 5vi as a function of the electron 
density for various potential amplitudes $* aui . The parameters are similar to those 
in Fig. 2. Insert: The measured velocity change of a SAW as a function of the gate 
voltage for different high-frequency (RF) powers applied to the IDT1; / = 114 MHz. 
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Fig. 5. The calculated intensities of higher harmonics with n = 2, 3 and 4 as 
functions of the potential amplitude $f aiu for two electron densities N s = 10 10 cm~ 2 
(upper part) and N s = 10 11 cm" 2 (lower part). The parameters are similar to those 
in Fig. 2. 

Fig. 6. The measured ratio T 1 P/I ae as a function of the high-frequency (RF) 
power P for / = 340 MHz; T = 300 K. The acousto-electric current I ae was 
measured in its maximum. 



27 



V 



contact 
IDT1 



ct 



v 

t 



v 




T 



contact 



IDT2 




semiconductor film with 
InCiaAs quantum well 



LiNbO, 



Fig. 1 



This figure "Fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0005 1 22 v 1 



This figure "Fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0005 1 22 v 1 



This figure "Fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0005 1 22 v 1 



1.0 




1.5 



in 

c 



Oh 0.5 P 



f = 340 MHz 
T = 300 K 



o r i 

50 100 150 200 

RF power P (mV) 



Fig. 6 



